Observation of breather-like states in a single Josephson cell 
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We present experimental observation of broken-symmetry states in a superconducting loop with 
three Josephson junctions. These states are generic for discrete breathers in Josephson ladders. The 
existence region of the breather-like states is found to be in good accordance with the theoretical 
expectations. We observed three different resonant states in the current-voltage characteristics of 
the broken-symmetry state, as predicted by theory. The experimental dependence of the resonances 
on the external magnetic field is studied in detail. 



I. INTRODUCTION 

Nonlinearity and discreteness in lattices lead to the 
existence of spatially localized stable solutions p], ^ [|. 
These solutions, called discrete breathers (DBs), are 
characterized by time periodicity and a strong spatial lo- 
calization of energy, on the scale of the lattice constant. 
The localized mode can be vibrational or rotational, de- 
pending on the lattice and type of excitation. In the sec- 
ond case the DB is called rotobreather. DB were found 
to exist in a large variety of nonlinear discrete lattices 
including weakly coupled optical wave guides Q , crystal 
lattice in solids 5], antiferromagnets Q and Josephson 
junction arrays 7, ||, 10, The latter system is 

a convenient experimental tool to investigate the exis- 
tence and dynamics of DBs. The discrete lattice here is 
constituted by an array of coupled underdamped small 
Josephson junctions (JJs). 

A JJ is composed by two superconducting layers sep- 
arated by a thin tunnel barrier. Its dynamics can be 
described in the resistively-capacitively-shunted-junction 
model (RCSJ) |Q. The equation for the superconduct- 
ing phase difference between the two superconductor lay- 
ers constituting the junction, ip, is written in the normal- 
ized form as 



7 = (p + atp + sin ip. 



(1) 



Here the time unit is normalized by the inverse of the 
Josephson plasma frequency lu p , 7 is the bias current / 
normalized by its critical value and a is the damping. 
This equation is formally equivalent to the equation of 
motion of a weakly damped pendulum in presence of an 
external torque 7. This mechanical analog allows for an 
easy description of the two possible states of the J J. If we 
apply a torque 7 > 1 the pendulum will start to rotate 
with an average frequency < ip > . After the pendulum is 
driven to the rotating state, it will stay in this state, even 
for 7 < 1, as long as the external torque is large enough 
to compensate the damping term. So that at some range 
of 7 both static and rotating states are possible, depend- 
ing on the initial conditions. By the second Josephson 
equation V = $o/2tt < (p >, these two states are the 
equivalent, respectively, to the superconducting and re- 
sistive states of the JJ. In experiment by measuring the 



voltage on the JJ it is possible to distinguish between 
these two states. 

Josephson ladders allow for the existence of both oscil- 
latory and rotational localized modes Jl3| . However only 
the latter ones are easy to detect experimentally. A roto- 
breather in the Josephson ladder is constituted by some 
of the JJs in the resistive (rotating) state, while the rest 
of JJs stays in the superconducting (oscillating) state. 
The rotobreathers were experimentally detected by mea- 
suring the I — V characteristics of the JJs in different 
points of the ladder |?], |l0j and also directly visualized 
by low temperature scanning laser microscopy j^, |^| . 

Josephson ladders can support linear cavity modes 
that can interact with the localized modes jlO], |ll|, |TJ, [l5| . 
The interaction of the linear modes of the ladder with the 
localized modes was found to lead to nonlinear features 
in the characteristic of DBs. Their complex dynamic in 
the region of the linear cavity resonances was investi- 
gated experimentally by the measurements of the I — V 
characteristics of various DBs . 

Single Josephson cells have been proposed Jl6|, |l7) to 
provide an insight into complex discrete breather states 
of larger systems such as Josephson ladders and two di- 
mensional Josephson junction arrays. Broken-symmetry 
states with characteristics similar to those of discrete 
breathers were proposed to exist (l^, |l7|, We will 

call in the following these broken-symmetry states as 
breather-like states. 

Superconducting loops with three Josephson junction 
have also been recently proposed for the realization of a 
qubit |l9[ |2fj|| . Such a loop has an important property re- 
quired for the qubit, namely two closely positioned min- 
ima in the dependence of its energy on the external mag- 
netic field |2lJ]. Note, however, that in this work we are 
interested into phase-rotational states, which are clearly 
different from oscillatory states discussed for qubits. 

In this paper, we study a single Josephson cell with 
three small underdamped Josephson junctions. Two 
junctions of same area placed along the direction parallel 
to the uniform bias current, Vi V2, are called vertical 
junctions. The third junction located in one of the two 
transverse branches of the cell, is called horizontal H, see 
Fig. I 

A parameter of cell anisotropy, rj, is defined as the ra- 
tio between the critical current of the horizontal junction 
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FIG. 1: a) Sketch of the single Josephson cell, with three JJs 
(indicated by crosses), two vertical Vi, V2 and one horizontal 
H. b) Optical image of our experimental cell with 4x4 fxm 2 
hole and three underdamped Nb/Al-A10a;/Nb JJs. 



and the one of the vertical junctions, rj = Ich/Icv- The 
parameter of anisotropy describes the coupling between 
the two vertical junctions, so that in the limit 77 — ^ 
the two vertical junctions are entirely decoupled, while 
the larger is rj the stronger is the coupling. The cell is 
called isotropic if all three junctions have same parame- 
ters, anisotropic otherwise. The system is uniformly bi- 
ased by two equal currents /, applied through the vertical 
branches, 1-2 and 3-4, Fig. [I]. As long as the uniform bias 
current / is smaller than the critical current lev, the 
junctions stay in the superconducting state, Fig. |^a. As 
the current / exceeds the critical value lev, the system 
switches to the homogeneous whirling state, constituted 
by both vertical junctions in the resistive state and the 
horizontal one in the superconducting state, Fig. ||b. 
Both these cases are symmetric states and there is no 
current on the horizontal junction, so that the two ver- 
tical junctions behave identically. We are interested in 
the last two cases, Fig. ||c and d, which arc the inhomo- 
geneous breather-like states. These two states have the 
same behavior and are characterized by the horizontal 
junction and one of the vertical junctions in the resistive 
state, while the other vertical junction remains in the su- 
perconducting state. In these cases, obviously, some part 
of the bias current flows through the transverse branch. 





M, |[. We bias the system by a local current, Ii oc , through 
the transverse branch 4-2, Fig. [j], and rise this current 
until the system switches to a breather-like state. After- 
wards, we increase the bias current / and simultaneously 
decrease Ii oc to zero, keeping the system in the generated 
state. Alternative procedure is first apply a homogeneous 
bias current / < lev at a value where we expect to find 
the breather-like state stable, then rise the local current 
Iioc to generate it, and finally reduce I/ oc back to zero. 
By using the two different directions of the initial local 
current it is possible to generate both breather-like states 
shown in Fig. g. 

Different states of the system can be experimentally 
detected by measuring the dc voltage on Josephson junc- 
tions. As for Josephson ladders, the broken-symmetry 
states are the consequence of the characteristic property 
of small underdamped Josephson junctions, that is the 
coexistence, for some ranges of the bias current, of two 
different stable states, superconducting and resistive. 



II. DYNAMICS OF THE SYSTEM 

The dynamics of the system are described by the phase 
differences across the three Josephson junctions, (pv l: 
ipv 2 ; fH ■ Using the current conservation law in each node 
of the cell and the flux quantization for the cell, one can 
get the equations of motion for the system jl7| f 



7 - irifVt - <Pv 2 + 2tt/) 

PL 

l+-S-(fv 1 -<PV 2 +VH + 27T/) (2) 
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FIG. 2: Possible states of the system, (a) superconducting 
state, (b) homogeneous whirling state, (c) (d) the two sym- 
metric breather-like states 



where the time is normalized by the inverse of the Joseph- 
son plasma frequency u> p . The influence of the external 
magnetic field is introduced by the term 27t///3l. The 
parameters of the system are the uniform bias current 
normalized by the critical current of the vertical junc- 
tions, 7 = 1 1 lev-, the normalized self-inductance of the 
cell, Pl = 2ttLIcv/$o (where L is its real inductance), 
the damping, a = \J <&o/(2ttIcR 2 C), and the frustration 
parameter / — <fr ext /$Q, that can be tuned by the ex- 
ternal magnetic field H . The operator M represents the 
current through a single junction in the RCS J model and 
is defined as: Af((p) = (p + aip + sin (p. 

Using a simplified dc model p7| , where the junctions 
in the resistive state are considered as resistances and 
the ones in the superconducting state just as shorts, one 
can estimate the dimensionless breather frequency f2 = < 
Lp >, that is in fact the voltage across a junction in the 
resistive state normalized by ^>olu p /2tt: 



The procedure to generate breather-like states in ex- 
periments is similar to that used for Josephson ladders 



a(l + 77)' 



(3) 
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From that one can also evaluate the region of existence 
of the breather-like state fl7| . The minimum current at 
which a broken-symmetry state may exist is 



4a . 

7r = —(! + //)• 

TT 



(4) 



At 7 < 7 r the system is retrapped to the superconducting 
state. The maximum current for a breather-like state is 



7c 



l + 2? ? 



(5) 



Above this current the system switches to the homoge- 
neous rotating state. 

The simplified dc model neglects the influence of the 
small oscillations of the Josephson phase that may cause 
an instability and nonlinearity in the characteristic of the 
breather-like state. Neglecting the damping and lineariz- 
ing the equation of motion, Eqs. (Q), around the breather- 
like state, Eq. (||), it is possible to calculate the two char- 
acteristic frequencies of these oscillations Jij], : 



where: F 
G 
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cos(ci). 



Here C\ is the phase shift of the vertical junction in 
the superconducting state evaluated by the dc model, 
ci = arcsin[7(l + 277) /(l + 77)]. The presence of the 
electromagnetic oscillations leads then to primary res- 
onances in the breather-like state, when the frequency 
of the breather VI locks to the frequency of the electro- 
magnetic oscillations oj + or and to parametric and 
combination resonances, respectively, at VL = nco±, where 
71 > 2 is an integer, and Q — n(cu + ± cj_), where n is an 
integer @ @. 

Using the rotational approximation j22j, two different 
dependencies of the resonant step amplitude A7 on the 
frustration have been obtained for primary and paramet- 
ric resonances Esl: 
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Expressions (^) and (||) have the same periodicity in / as 
the critical current in the whirling state, but, depending 
on the parameters of the system, they may be asymmetric 



with respect to / = 0. Moreover, the external magnetic 
field can considerably change the amplitude of the reso- 
nance and, in the case of primary resonance, even make 
it to vanish for some values of the frustration. 



III. EXPERIMENTAL RESULTS 

Our experimental system is a single cell with the hole 
area 4x4 /1m 2 , containing three small underdamped 
Nb/Al-AlO^/Nb Josephson junctions ||, see Fig. g. 
The critical current density, measured at 4.2 K, is Jc = 
133 A/cm 2 . The area of each vertical junction is about 
7x7 /1m 2 , while the horizontal one is 3 x 4 /1m 2 . The 
anisotropy of the system, evaluated from the design, is 
then 77 = 0.24. The parameter of self-inductance of the 
cell at 4.2 K, evaluated by measurement of a SQUID of 
same geometry, is (3l — 1-47. The two equal bias currents 
are introduced in parallel using two large external resis- 
tances of 1.5 kfi, so that the bias current remains uniform 
also in the broken-symmetry states. The damping of the 
system at 4.2 K is evaluated from the retrapping current 
of the homogeneous whirling state, a ~ (ir/A)(I r /2Icv), 
that yields a ~ 0.03. The measurements are performed 
(6) in a 4 He cryostat [24j in the temperature range between 
4 K and 8 K, so that the parameters of the system Jc, 
a, Pl and u p can be varied. For each temperature the 
damping a and the critical current density Jc are evalu- 
ated from the characteristic of the homogeneous whirling 
state, while the parameter of self inductance /3l is evalu- 
ated from the critical current density. Finally, the value 
of the plasma frequency is obtained at every temperature 
from fitting the voltage of the parametric resonant step 
to the expected value 2uj + . 

For convenience of comparing with theory, the experi- 
mental data are shown in the previously introduced nor- 
malized units of the bias current 7 = I /lev, breather fre- 
quency fl = V/ (uj p &o/2ir) and frustration / = & ext /$o. 

In the range of temperatures from 4.2 K to 7.75 K 
(0.02 < a < 0.35) it was possible to excite both breather- 
like states and they were found stable in a finite range of 
the bias current. No breather- like state was found stable 
for temperatures above 7.75 K (a > 0.35). 

The minimum bias current below which the breather is 
retrapped increases linearly with the damping, see Fig. ||, 
in good agreement with the theoretical expectation of 
Eq (^). The maximum current that allows for the ex- 
istence of the breather-like state slightly decreases with 
increasing the damping, see Fig. ||, whereas from Eq. (|^) 
it is expected to be constant (j c = 0.84). 

In fact, the temperature variation leads to an alter- 
ation of several parameters, see Fig. |4| In particular the 
decrease of the self inductance Pl at high temperatures 
leads to stronger interactions of the breather-like state 
with the small oscillations. This effect may strongly af- 
fect (in the range 0.1 < (3l < 1) the region of existence of 
the breather-like, especially its upper border fl7|| . In fact, 
the experimental data show deviation from the expected 
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FIG. 3: Maximum (solid squares) and minimum (solid dia- 
monds) currents for the breather-like state. The solid lines 
show the theoretical expectation of the dc model, Eqs. (^) 
and (0). 



a 0.2 




FIG. 4: Dependence of the damping, evaluated by the retrap- 
ping current (solid squares) and by the sub-gap resistance 
(solid circles), and of the self-inductance (open triangles) on 
the temperature of the system. 



value for a > 1.15, which is in the range of the self induc- 
tance between 0.9 and 0.3. No resonances were observed 
in the I — V characteristics at temperatures below 6 K. 
The small damping of the junctions at this temperatures 
seems to be unfavorable for reaching the small frequencies 
of the electromagnetic oscillations where the resonances 
are expected (fl m in — 5.5). Increasing the temperature 
of the system had two helpful effects for seeing the res- 
onances, namely the increase of damping and decrease 
of self inductance, see Fig. [|. As the consequence of 
the larger damping, the breather-like states were stable 
down to smaller frequencies, whereas a smaller induc- 
tance moved the frequencies of the resonances up. 

At the temperature of T = 6 K it was possible to 
observe the first resonant step. The parameters of the 
system at this temperature are a ~ 0.11, J c — 92 A/cm 2 , 



(3l = 1-02 and uj p /2n = 37 GHz. This resonance is a step 
in the characteristic, shown in Fig. |^. Comparing this 
resonance with that at higher temperatures and with the 
theoretical expectation we deduce it to be the parametric 
resonance at 2w+. 
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FIG. 5: Measured current-voltage characteristic of the 
breather-like state at damping a ~ 0.11 and frustration 
/ = —0.1, where the magnitude of the resonant step is maxi- 
mum. The dotted lines show the expected frequencies of the 
electromagnetic oscillations, Eq. (^). 

At the temperature of T — 6.65 K we detect two addi- 
tional resonances in the characteristic of the breather- like 
state, see Fig. ^. At this temperature the parameters of 
the system are a ~ 0.18, J c = 70 A/cm 2 , (3l — 0.77 and 
(jjp/2ir = 30 GHz. The two larger resonances, A and B, 
are different by the integer factor of two, as expected for 
2lu + and u> + . 

The dependencies of the current amplitude on the ex- 
ternal magnetic field for steps A and B are strongly differ- 
ent. The resonance A has its maximum amplitude close 
to zero frustration, see Fig. |^a, and disappears in a wide 
range around / = 0.5, see Fig. (?]. An other resonance B 
is always present but has nearly opposite behavior. At 
the top of the resonance B, around frustration / = 0.4 
or / = —0.6, the system shows an instability by switch- 
ing either to the homogeneous whirling state or to the 
superconducting state, as shows Fig. |7|. In this range 
of frustration the measurements of the resonant step B 
were performed with less accuracy and faster because of 
its poor stability. We argue the larger step B is due to 
the primary resonance at £1 = lu + , whereas the step A is 
due to the parametric resonance at f2 = 2lu + . The third 
smaller nonlinearity, C, is at a frequency between the 
other two, close to the theoretically expected value for 
the composed resonance atO = w + +w_, see Fig [| This 
resonance was only present in a small range of frustra- 
tion around zero, where from the theory, the composed 
resonance is expected to be maximum. Because of its re- 
duced amplitude this resonant step has no hysteresis and 
its dependence on the magnetic field was not measured. 
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FIG. 6: Measured current-voltage characteristic of the 
breather-like state at the damping a ~ 0.18 and for:(a) 
/ = —0.1, (b) / = 0. The dotted lines show the expected 
frequencies of the electromagnetic oscillations, Eq. (^). The 
insets show a zoom into the region of interest. 



FIG. 7: Measured current-voltage characteristic of the 
breather-like state in presence of instabilities at damping 
a ~ 0.18 (a) / = -0.6, (b) / = 0.4. The dotted lines show 
the frequencies of the electromagnetic oscillations evaluated 
from Eq. (§). 



In Fig. gp, c we present measurements of the current 
amplitude of the resonant steps A and B on the exter- 
nal magnetic field. The parametric step at T = 6 K, 
see Fig. [| is higher than that at T = 6.65 K, Fig. ||. 
Moreover, because of the larger damping in the latter 
case the step vanishes in a larger range of the frustration 
(-0.95 < / < -0.25 and 0.05 < / < 0.75). Therefore, 
in order to observe a more clear modulation of this reso- 
nance, we report in Fig. ||c its dependence on the external 
magnetic field at T = 6 K. 

As expected from the theory, Eq. ((?]) and Eq. (||), the 
magnetic field dependence of the resonance amplitudes 
show an asymmetry with respect to / = 0. Note, that 
such asymmetry is not observed for Ic(f), see Fig. ||(a). 
The asymmetry depends on the parameter of the system 
and particularly on the bias current where the resonance 
occurs. This behavior can be qualitatively attributed to 
the current flowing in the transverse branch of the cell 
in the broken-symmetry state. The experimental data 
are in rather good agreement with the theoretical expec- 



tations. We fitted the maximum and minimum of the 
resonant steps in order to find the proportionality coef- 
ficient in Eqs. (@) and (§). The small variations of the 
experimental data around the maximum of the primary 
resonant step could be due to the appearance of insta- 
bilities, as discussed in relation to Fig. |?]. This behavior 
was also found in simulations [ fl8| . In the case of para- 
metric resonance the presence of a range of frustration, 
that increases with the damping a, where resonant step 
disappears, is also expected Q. Note, that the region 
where the experimental data show the largest deviation 
from the expected behavior is nearly the same for the 
two cases (0.3 < / < 0.6). 



IV. CONCLUSIONS 

We observed breather-like states of broken symmetry 
in an anisotropic single-cell system. For temperatures 
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lower than 6 K the current- voltage characteristics of the 
breather-like states were stable and did not show any 
resonances. Increasing the dissipation by rising the tem- 
perature of the system we decreased the frequencies of 
rotation of the breather-like state in such a way that it 
matched the frequencies of the resonances in the cell. In 
this region, the primary, composed and parametric reso- 
nances were observed. As expected, in all three cases, the 
interaction of the electromagnetic oscillations with the 
breather-like state leads to resonant steps in the current- 
voltage characteristic. The large steps due to the primary 
and parametric resonances show hysteresis, whereas the 
composed one induces just a small nonlincarity in the 
characteristic. The measured dependence on the external 
magnetic field of the magnitude of the primary and para- 
metric resonant steps is in good agreement with theory. 
In the case of primary resonance, the external frustration 
can strongly increase the amplitude of the resonant step 
and lead to instability of the breather-like state. The 
composed resonance, for this values of the damping, was 
present just in a small range around zero frustration and 
did not lead to instabilities. The region of existence of the 
breather-like state was also measured and found in good 
accordance with the expectations. At smaller values of 
the self inductance parameter, the stronger interaction of 
the breather with the electromagnetic oscillations leads 
to a contraction of the region of existence. 




FIG. 8: (a) Dependence of the critical current of the homo- 
geneous state on the frustration, (b) , (c) Dependence of the 
magnitude of the resonant step on the frustration for primary 
resonance A (measured at T = 6.65 K) and parametric res- 
onance B (measured at T = 6 K) respectively. The solid 
lines are the theoretical curves given by Eq. ([?]) and Eq. (j§|), 
respectively. 
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